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ABSTRACT 


Modern naval vessels, which depend more on a reduction in 
the displacement per shaft horsepower than an increase in the 
Peeeyarag ratio to improved performance, must increasingly 
turn to propulsion systems with characteristically high power 
densities. Systems planned for or currently employed in high 
performance craft all have inherent disadvantages which 
reguire new approaches. One such approach is the use of 
electrical propulsion systems which allow operational flex- 
ibility and an ease of power transfer not found in many other 
Systems. The utilization of superconducting materials in the 
field windings of propulsion machines significantly improves 
the power density of the system while retaining other advan- 
mereeous Characteristics of electrical propulsion systems. 


Many superconducting motors either in the conceptual or 
prototype stages of development, require the superconducting 
Field winding to transfer machine torgue to the naval craft 
Pmeeley 1 turn, requires a Substantial structural interface 
between the superconducting element and the environment. 
such a structural interface results in a significant heat leak 
into the cryogenic components reducing system efficiency. 
Other designs place the cryogenic element on the rotor, 
h-eurcing a rotating liquid transfer coupling and all of the 
technical complexities that such an element entails. 


This thesis proposes the use of a dual field motor; the 
main power producing field to be superconducting, and an ambient 
temperature concentric compensating field configured to remove 
the load torque from the cryogenic element. Both fields are 
WOUMomenetne Stator, thus eliminating the liquid transfer pro- 
blem; the three-phase armature 1S wound on the rotor. A 
mathematical model is developed and the size and parameters 
for a 22,700 horsepower motor are determined. The resulting 
model is analyzed and various modes of control discussed. 
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i LNTRODUCTION 


tiemiese On Ssujerconducting field windings in synchronous 
machines iS now an accepted and experimentally substantiated 
technique, and it is evident that these machines have the 
greatest potential applications where weight and power density 
are Of significance in the system design philosophy. Because 
of the inherent weight and volume savings and system flexibility 
associated with superconducting power generation and 
utilization, superconducting machines present an attractive 
alternative in the design of naval propulsion systems. 

As progress is made from lower to higher performance 
water-borne vehicles, increases in Speed are entirely attri- 
butable to increases in the shaft horsepower per ton displace- 
ment rather than to improvements in the vehicle lift-drag 
merase in addition, hull structure and propulsion plant 
together have in past designs accounted for an average of 60% 
of the light ship weight, and future marine engineering 
studies, especially those in the naval systems domain, should 
concentrate on such developments as superconducting machines 
which promise significant weight and volume Savings and 
superior acoustical performance. Electric power trains, 
beemise ot their flexibility, are particularily attractive in 
advanced ship concepts such as SWATH, hydrofoils, and surface 
effect vehicles which require power transmission in strictly 
eenmelnec scaces and Lransmission trains which change direction 


abruptiy. 





All of the superconducting synchronous motors now either 
in the conceptual or prototype phases of development, utilize 
a rotating superconducting field winding which requires 
@evogenic fluid transfer to the rotating member. Many designs 
also require the Superconducting field winding to transmit 
machine torque to the ship hull which requires a substantial 
structural interface between the cryogenic elements and the 
Ship's hull. This results in significant thermal leakage to 
the low temperature machine components. 

Dieemeter evaluated herein will have a stationary super- 
conducting field winding and an ambient temperature rotating 
armature winding, the armature being the central, interior 
member (see Figure 1). The field winding iS Surrounded by a 
Solid, iron flux return path, operating in a magnetically 
saturated condition. A third direct current ambient tempera- 
ture winding is provided to assume the load torque and 
relieve the superconducting field winding of this duty. This 
design eliminates the need for a rotating cryogenic transfer 
Peupling and will substantially reduce the structural inter- 
face between the superconducting field and its environment. 

The compensating field winding, if feasible, may be 
located either between the armature and main field or between 
the main field and shield. This thesis will examine both 
configurations and draw the appropriate conclusions from the 


Mmiomveteakt results. 
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This paper developes the motor in a conceptual framework 
with emphasis on naval applications and high performance 
vehicles. Figures 2, 3, and 4 depict three possible applica- 
tions to naval ship systems that are currently considered 
"state-of-the-art" to the naval architect and marine engineer. 
The motor developed herein may find useful applications both 
in these systems and many of those to be conceived in the 


fucure. 
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ieee HYSICAL DESCRIPTION OF MOTOR 


Figures 1 and 5 depict the approximate physical con- 
figuration of the compensated motor. The illustrations used 
throughout this paper address only the two pole-pair (four 
pole) machine, but the analysis herein applies to configura- 
tions with any desired number of pole pairs. 

The three-phase armature ee ound On “che sO LOr land 
supplied power through a high voltage slip ring assembly. 

The armature is to operate at ambient temperatures and is 
cooled by conventional non-cryogenic techniques. It is 
assumed that the upper constraint On armature current density 
[/emapproximately 2.5 x 10° amps per square meter. 

The superconducting field winding is positioned either 
outside of the compensator winding or between the armature and 
compensator windings. It 1s to operate at a temperature less 
than 5°C. As discussed earlier, the field winding may be 
fixed to the machine housing through an insulating medium 
Wamem also functions as the supporting structure for the 
field. The field structure iS assumed to contain little, if 
any, ferromagnetic material. 

An idealized depiction Steeve fFiela winding Support 
meeterice 1S Shown in Figure 6. This configuration is 
particularly applicable to a machine with the compensator 
within the field winding. Oe neces the field could be 
peperted radially and axially through the end housings of 
miewmachine. The Optimum support structure could be depicted 
as some configuration of bicycle spoke supports made from an 


Seeective thermal insulator. 


23 





"ANV1d Z-A‘YMYOLOW ONILONGNOOMZdNS 
Q3LVSNSdWOO-ANOYOL 4JO NOILOSS SSOHD °S 3uNdIS 


MISWASSV 


ONIN diqs 
STWNIWYSL YSMOd 


ONY SINSIOAND i “em a a 


LSVHS Sg et 
yaiadOud 
AISWASSY 


SYNLVWYYV ¢ N ONIEY Ag 


SNIONIM 
Olalsa SNILONGNOOYAdNS 


ONIGNIM Oe ae 
HOLVSNSdWOOS LNZISWY : 
ONISNOH YOLOIN we 


ONV QISIKS SILANSVW 





24 





"INSWSAITS Aisid NIVW 
ONILONGNOOYSdUNS YOS SYNLONYLS LkOddNS G3ZITVaGI “9 SYNdls 


SJHNALONYLS 
LYOddNS G3LVINSNI 


Gilaiizs 
Q ald 






uy) 
on 








iiewecompensating Winding 1s physically located either 
outside of the field winding or between the field and the 
armature windings. It 1S to operate at ambient temperatures, 
meme COOled by non-cryogenic techniques, and to contain 
meeele, if any, ferromagnetic material. The current density 
is constrained at a maximum value of 2.5 x 10°° ampS per 
eoeatre meter. | 

The motor examined herein iS envisioned to be designed 
as a modular element in a propulsion system which must 
operate with a minimum down time after failure. 

The motor will be fitted with easily disconnected 
meeeings in Order to facilitate rapid replacement upon failure 
at shore based facility. The cryogenic auxiliaries as shown 
in Figures 2, 3, and 4 are also to be modular packages which 


may be rapidly replaced by new or refurbished units. 
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meee GORE TTCAL CHARACTERISTICS AND RATING EXPRESSIONS 


An examination of Figure 7 reveals the requirement for 
zero torque on the superconducting field winding in terms 
of the machine's magnetic fields. The fields are referenced 
to the stator, and the resulting field vector at any location 


may be found as: 


pees B, + B. (1) 


where the resultant of the compensator and superconducting 
field winding fields is expressed as an effective field flux 


density: 


ne eeond Bo are the fielda vectors of the armature, 


ey F 
compensator and field respectively. 

A requirement for zero torque on the superconducting 
meld winding is that the vectors Be and Be be colinear as 
reflected in Figures 7 and 8. 

The following constraints placed on machine design 
parameters are specified as: 

l. Maximum field current limited to 1.25 x 10°%a/m* 


2. Maximum armature current limited to 2.5 x 10°%a/m? 


Sa coo mum Gompensator current limited to 2.5 x Oley ie 


Z| 
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4. Machine magnetic shear stress is to be 
consistent with current design practices 
Reference 6 presents a detailed discussion of constraints 
1 through 3 above; and reference 7, Figure ]1, shows magnetic 
shear stresses found in recent designs. 

Referring to Figures l, 7, and 8, the angle a represents 
the angle between the axis of the Superconducting field and 
the compensator field.Figures 7 and 8 show all possible 
operating loci for the machine field vectors constrained as 
above and operating as a motor. If angle a is decreased 
from 90° electrically to some other fixed position, the 
domain of operation decreases also. A comparison of Figures 
7 and 8 shows this clearly. If a is permitted to vary 
rather than fixing the relative position of the two fields, 
an additional degree of freedoin is introduced into the 
analysis. Figures 9 and El represent the resulting voltage- 
current relationships for the compensated machine. AS a 
result of the colinearity requirement imposed on the fields 
ae and De, E, and V,, must pib-omibewcolinear as. shown in 
Figures 9 and El. This results in a somewhat simplified 
analysis as presented in Appendix E. The analysis in 
Appendix E treats a aS a constant < 90° electrically. 

Meterming to Figures 9 Or El, the internal excitation 


voltage arising from the field vector B. is expressed as: 


im 


= _ ol 
Bi = GU era cos W cot a) (ae) 
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FIGURE 9. VOLTAGE-CURRENT RELATIONSHIPS FOR OVEREXCITED 
AND UNDEREXGITED GOMPENSATED MOTORS. 
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Mepbizing this relationship and the appropriate expressions 
her Machine inductances, machine power in volt-amps is 


found to be (see Appendix FE): 








0 ale V 
24 Dawe \ PP ewa, ,_ Powe Jb 
Pe= —— WwW UdSJS-2 (l-y ie CG 7S Ue Sint ) (=) 
fan & ° a i rn EO. =LOE 2 io 
(E8) 
when 
Ver 
—— = 1-x (sin W - cos wp cot a) (E14) 
a 


and, Xt the internally normalized synchronous reactance is 
given by equation E16. Equation C38, the armature torque 
expression, may be written as 
= A_[{A.si c + 
Es if 5 Sin k @ A 


3 J sin k (d-a@) ] (E20) 


where Avs Aas and A. are defined in equations E17 through 

moe and are functions of Jos Jer and machine geometry. 
The armature torque, Toe may also be expressed by 

Pem@erton b2l and the resolution of equations E20 and E21 


results in a quadratic expression for the angle between the 


Sebo, Ot phase a and the superconducting field (angle 9): 
in? (k in(k = 
On sin* (ko) + om Sin(kdo) + oe 0 22) 


where 0., QO and On ane Getined in equations E23 through E25; 


Ee 


and are functions of Je a Jae and machine geometry. 


ie 


32 





Utilizing the constraint expression for zero torque on 
mene superconducting field winding given by equations C38, 


moor 26: 
Sesi(k¢d) = Je A cos (ka) (226) 


A fourth order equation for the compensator current density 


is found in Appendix E to be 


et NJ + A 


2 = 
oe 4 346 1° A,J . te Ay = 0 (E27) 


where Ay Eineoue lh A. 


The coefficients are also given in equation E34 for the 


are defined in equations E29 through E33. 


specialized mere with four poles (p=2). 

Using the above expressions in light of the aforementioned 
constraints, a FORTRAN program, THESIS I, was developed as a 
design tool to facilitate the selection of a feasible machine. 
"Feasible" is defined here to indicate a machine that meets 
the size, shear stress, and current density constraints 
discussed earlier. Other machine operating characteristics 
Buen as per-unit synchronous reactance, terminal or phase 
voltage, line current, and power factor are also considered. 
Appendix F presents more detailed information for use of the 
meOogram. The program Biaees Bn incremental adjustment of 
the Bete o and determines the required compensator current 


Gensity to maintain zero torque or the superconducting field. 
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As indicated in Appendix E it is also possible to approach 
the problem by incrementing J and determining the angle a 
required for zero torque (see equation E47). 

Because the solution of equation E27 yields four roots, 
some complex, some real, the valid real root may be found by 
Mi@iestituting back into equation E22 and then E26. 

Bagures 7 and 8 indicate Bree anclesa Sel at 90" 
electrically allows the greatest possible range of machine 
operation. Of course, by setting a to any value requires an 
adjustment in the compensator current density in order to 
maintain a net zero torque on the superconducting field. fThis 
particular mode of operation permits the designer to use 
fixed, rigid interfaces between the compensator and super- 
conducting field winding and the supporting environment. 

An attractive alternative to this technique is to permit 
the angle a to seek its own natural value as power level 
changes and J. is held constant. This requires a flexible 
interface for cryogenics and electrical power cables but 
eliminates the need for an external control loop. A third 
possible, but less attractive method, is to fix the compen- 
sator current Jas and mechanically position the field or 
compensator through a feedback control loop. 


Equation E40 
p= KJ - K J, ° (sin Pees UecOt 0.) (E40) 


K, and Ky defined by equations E36, E38, and E41, indicates 
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me@emutility that angle o may offer in adjusting the machine 
ewer level. It is noted that if ao is used to adjust power 
level, J Moco vil rca tO Malntain the zero torque condition. 


TE 1S also possible to express the field current density, 





Jer as: 

J = K,. =f K,d_ (sin eqs y GC 30 ) (E42) 
where K. 1s defined by equation £43. If A is defined as 

m= sin WW - cos W cot a, (3) 
E42 may be written as: 

— ey) Je = ep (B44) 


Figure 10 depicts A as a function of the power factor 
angle ~ for various values of a. If the machine is over- 
excited (quadrants I and IV of Figure 10), » may have either 
positive or negative values, depending upon the value of 
angle a. If the machine is underexcited (avadrants II and III 
of Figure 10), A may only assume negative values. Figure ll 
depicts plots of equation E44 for a fixed power factor angle 
with the machine underexcited. An examination of Figures 10 
and 11 shows that the slope of equation E44 may have either 
positive or negative values depending on the angle a for the 


overexcited machine. Conversely, Figures 10 and 12 indicate 
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only positive slopes for equation F44 if the machine is under- 
feeertcd. A classical "V-curve’ for the synchronous motor can 
only be obtained if a = 90°. From Figure 10 this yields the 
same magnitude of A regardless of the machine excitation, 
poea “V-Curve" is obtained for each value of +tW. 

Reeerwise, 1f a 7 90°, frompPigure 10, it is evident 
that the magnitude of A is not the same for both a + w and 
eC Therefore the "V-curve" will not be symmetric about 


peevertical Jine through J 5 = Kee 


Sey, 





ive RESULTS 


erbiging tae program Listed in Appendix E (THESIS I) a 
feasible machine, whose physical and operating parameters are 
listed in Table 3, was developed. This machine employs a 
compensator physically positioned concentric with and outside 
of the main superconducting wae winding. 

Figure 13 is a plot of some of the designs investigated 
during the course of the development of this thesis. All 
machines in Figure 13 have the compensator winding physically 
located between the armature and the main field windings. 
Adjacent to each data point is printed the per unit synchronous 
reactance at the power level indicated. Each machine plotted 
in both Figures 12 and 13 are designed at the upper limit of 
armature current density (gi = 2 Moceaym- wk) machines 
plotted in Figure 13 were designed at the upper limit of field 
Seeerent (1.25 x 10°a/m*7). All machines were assumed to 
operate from a 60hz bus and configured for two pole-pairs. 

Piamce sl > ts a Plot Of magnetic shear stress as a 
function of power for a number of machines with the compensator 
Smaesidqe of the main field winding. 

Tables 1 and 2 are provided to augment the data plotted 
Pmeagures 12, 13, and 14. 

The machine presented herein is rated at 22,700 horse- 
power before losses, and is intended to operate as part of 
the main propulsion system for the aforementioned hydrofoil. 


ii@esenotce 15 purely arbitrary, and the SES or SWATH 
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meeerecation could have been chosen for illustration with 
correspondingly higher horsepower ratings. 

Figure 16 shows Jo go 2 funce lon of J fOmea faxed 
field current and a fixed value of the phase angle between 
the main superconducting field and the compensator field. 

@Mhne plots indicate a relatively simple control problem, with- 
Out the need for mechanical components in the control loop. 

As discussed earlier, another option available to control 
the torque on the main field winding is physical angular 
adjustment of the main field winding (angle a in Figure 1), 
with a constant value of the compensator current density, un 
Pigure 17 shows the effects of varying angle a and J. on 
machine power output. The small circles on the a = 10° and 
mee splots indicate the operating point at which J. exceeds 
Biemeonsctraint of 2.5 x 10°a/m*. In the other plots, a 2 30°, 
mies constraint on J. is not exceeded, but as shown,a similar 
Bomstraint on oa 1s reached prior to reaching the we GOnSseraiIn. 
Therefore, for the machine described in Table 3 the limiting 
value for the current density in the compensator is not exceeded 
mora 2 30°. 

Figure 18 depicts the variation required in angle a for 
Bemstant values of J. as J, Or power varies. It is again noted 
that the power level at any angle o and any armature current 
density, Jot He olcOeamLiumcelOnJOL sLle compensator current 
density, Jy. In the limiting case of a = 0° electrically, 
the compensator would act as only as a power augmenting field 


element and in order to maintain zero torque on the 
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3.2 RATED POWER 
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Mire mieqmnceld, en infinite amount of current would be 


xed in the compensator. 








a 
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PeeeeDLOoCUSSION OF RESULTS AND CONCLUSIONS 


AS indicated in Figures 13 and 14, a significant number 
of machine Aste iee investigated uSing the analysis 
developed and constraints specified herein. THESIS I, as 
listed in Appendix F, was employed extensively as a design 
EOoOl. 

No machine with the compensator located between the 
superconducting field and armature was found to be satisfactory 
Meee because of a firm constraint violation (high current 
densities or low shear stresses) a because of unacceptably 
low values of per-unit synchronous reactance. It is reco- 
gnized that low values of per-unit synchronous reactance 
present problems primarily during machine transients and 
because of the extraordinary sequence of events that must 
occur to induce a major caSuality to the motor, low values 
may not be meaningful either to the designer or operator. 
Figure 13, auaqmented by Table 1, does indicate some machines 
feuch as machine A, D, E, or F) which are theoretically 
plausable but have either low per-unit synchronous reactances 
Or poor power factors at designed power levels, or both. The 
inability to obtain satisfactory machines at high power levels 
(2 20,000 hp) is attributed to the pronounced decrease in 
mierda strength in the radial direction. Consequently, with 
the compensator located between the superconducting field and 
the armature, lower power must be expected than that produced 


from a machine with the compensator located outside of the 


Del 





@reld winding, for the same superconductor, compensator, 
and armature winding areas. 

In summary, the machine with the compensator located 
inside of the superconducting field presents many design 
problems because each machine investigated at rated power 
either exceeds or approaches the design constraints. 

Smecne Other hand, a motor with the compensator located 
between the superconducting field winding and the shield is 
less demanding in terms of the deSign constraints. Figure 14 
depicts rated horsepower as a function of compensator current 
density for a number of designs investigated. The machine 
parameters used as constraints are not approached, and as 
mmerecated in Figures 14 and 15 and Table 2 acceptable values 
of per-unit synchronous reactance and power factor are easily 
realized. If Figures 15 and Dl are eelipared, Ene values of 
O for those machines investigated are acceptable. 

Table 3 presents data for a feasible machine to be 
Brolioved as a Main propulsor for a hydrofoil. The compensator, 
Superconducting field, and armature current densities are all 
well within the constraints outlined in Section III of this 
paper. An optimun machine was not investigated, and possibly 
a machine of smaller dimensions could be synthesized. 

Of the three options discussed earlier to control the 
torque on the Superconductors, fixing the angle a at 90° 
electrically PopeeOneenen autho Ss POInt Of View, the better 
choice because such a mode of operation allows the widest 


possible locus of operation as depicted in Figure 7. 
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Meee 3S. DATA FOR A FEASIBLE TORQUE-COMPENSATED MACHINE 


Parareter 

Design Horsepower 

Reyer Factor 

KVA Rating 

Design RPM 

Rated Compensator Current Density, J 


Rated Field Current Density, Je 


Rated Armature Current Density, J 


Rated Terminal Phase Voltage 


Rated Terminal Line Voltage, Vin 


Phase Angle Between Main and Compensator 
Fields, a 


Armature Winding Angle, oF 3 


Field Winding Angle, 6,, 


Compensator Winding Angle, Oc 


Armature Space Factor, A 


Armature Winding Inside Radius, Ro; 


Armature Winding Outside Radius, Roo 


Field Winding Inside Radius, Res 


Field Winding Outside Radius, Rey 


Compensator Winding Inside Radius, a 


Compensator Winding Outside Radius, Ee 


Shield Inside Radius, ee 
Shield Outside Radius, ee 


Shield Permeability, te 
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mebbs 3 continued 
Poaranecer 
Total Machine Length, a. 
Berpaight Section Length, & 
Synchronous Reactance (per unit), Xg 
Frmature Conductivity 
End Turn Length (Ro; qe Ee 
Effective Lengths For: 
self Inductance, Armature 
self Inductance, Field 
Self Inductance, Compensator 
fatual Inductance, Field-to-Armature 
Mutual Inductance, Compensator-to-Armature 
Conductor Loss 
Bady-Current Loss 
Losses at Rated Conditions: 
Conductions (compensator and shield) 
Eody-currents 


Mortal, percent of rating 


Ratio of & /D (see Table D1) 
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Bagure 8 indicates that for any other fixed value of a, a more 
limited locus 1S to be expected. Figure 16 shows that as 
angle ® is reduced from 45° mechanically for a four pole 
machine, a higher value of compensator current is required 


Memeny Particular value of J, (or any particular power level 


A 


since Vin is assumed constant). Such an arrangement also 
eliminates the need for moving interfaces between the sSuper~ 
conducting field winding or the compensating winding. This 
means that rotating fluid cryogenic couplings and slip rings 
are not necesSSary; and conseguently should both improve 
machine reliability and simplify the mechanical design. 

The primary disadvantages of employing this technique 
is a slower compensator response to machine load changes and 
the added acmauenikee of an external feedback loop required 
fee Control es 

Figure 17, plots of motor horsepower as a function of 
J, for the machine cesScribed in Table 3, indicates an 
increase in horsepower due to the improved armature-to- 
Sempensator coupling as angle © is varied. As discussed 
Garlier, if angle © is permitted to go to 0°, the compensator 
field aligns with the superconducting pmeildsand- acts Only as 
a power element, no longer capable of compensating the super- 
conducting field without infinite currents in the compensating 
winding. Some estan OF Course, would be realized at 
a = 0° betes the compensator would absorb some torque, but 
its value would be insignificant compared to that carried by 


the superconducting field. 
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A viable alternative, as briefly addressed above, but 
not as advantageous from the standpoint of hardware reliability 
and mechanical design is to mount either the superconducting 
field winding or the compensating winding on bearings 
@eestoily a self-aligning ball or cylindrical bearing) and 
merow the element to seek the proper value of angle a at some 
fixed ay as power level changes. Figure 18 predicts the 
behavior of the element for various values of J as the power 
level varies. The machine depicted is that described in 
Table 3, and the terminal voltage is assumed constant. The 
advantages of such a mode of operation are rapid response of 
the compensator to changes in load, and the elimination of 
the external feedback control system necessitated in the 
system Meo iloned above. 

Since either the superconducting field winding or the 
compensating winding can be gimbled, the designer has an 
Meeeeronal choice of configurations; both of which require 
moving interfaces with the environment. If the superconducting 
Meera 16 fixed, the rotating cryogenic coupling is not 
required and only power slip rings are necessary to provide 
power to the compensator which, of course, would be gimbledad. 
Otherwise, if the compensator is fixed, the main super- 
conducting field requires both cryogenic and electrical 
interfaces. Because of eee conclusion presented 
herein, en placed the compensatOr outside of the super- 
conducting field, a fixed compensator and a gimbled main 


field would further complicate the machine's mechanical design. 


56 





As may be seen above strong cases for either mode of 
Operation can be made; but because this machine is envisioned 
to be employed in naval or possibly marine vehicles, the 
author feels that machine reliability should play an 
Piportant role in the design decision. Therefore, a fixed 
Compensator and superconducting field with a variable 


compensator current 1S recommended. 
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VI. RECOMMENDATIONS 


Additional areas that could be investigated using the 


material developed in this thesis are: 


ae 


k= 
° 


Develop field equations in terms of Oo and 1 
and apply results to develop steady state and 
transient machine performance. 

Senseruce aS an initial prototype, a small non- 
cryogenic machine and examine machine behavior. 
On a systems level, develop alternative power 
trains based on system volume, weight, arrangements, 
efficiency, and complexity for various naval 

op reat LOns. 

Synthesize a control system model for the motor 
operating at a constant value of a, varying ae 
with power level. 

Develop transient and steady-state behavior of 


machine with variable shield permeability by 


modeling the shield magnetization characteristics. 
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APPENDIX A 


FIELD ANALYSIS 


This Appendix addresses the derivation of the basic 
field equations for the machine under consideration. The 
technique employed is that discussed in references 2 and 3 
which utilizes axially constant boundary conditions expressed 
as ratios of radial and azimuthal fields. The use of these 
ratios, referred to as surface coefficients, simplifies the 
arithmetic involved and permits a wide range of boundary 
conditions to be applied to the machine problem after the 
solutions are formalized. 

In the most general form, Ampere's circutial law may 


be expressed as: 


| 
~* 
I 
II 
cal 
o- 

ad ot 
> 


The primary concern here is with the fields within the 
machine air-gaps which are current-free, where Ampere's Law 


may be expressed as follows: 


=- = aD 
= A2 
V x H E (A2) 
Assuming that the E fields do not vary with time, 
Ampere's Law becomes: 
vaexH = 0 (A3) 


Se, 





The electromagnetic fields may be specified by a system 
@Sreequipotentials which are cut at right angles by lines of 
magnetic force. This property of orthogonality between the 


lines of force and the egquipotential surfaces leads to some 


relatively simple geometric conditions. Define: 
te magnetic flux function 
@? = magnetic potential 
The family of curves, % = constant and wW = constant, 


are orthogonal and this constraint may be expressed in polar 


meeroainates as: 


od? op _ _ d% dy 
dx ° Or ~ ~ rod roo 
From the theory of functions of a complex pedeiseee! 
the functions ¢ and W are analytically related via the 
complex variable 
Z=x + jy (A5) 


Genverting AS into polar coordinates 


m= COSG and y =r sing 


therefore 
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ot ay = £ir(coseé + jsiné)] (A6) 


Potterentiating AG partially with respect to r and 6 


successively, a comparison of the results indicates 








0 ,) 

ee -. (A7) 

0 > ,) 

ee > sr oo 
Differentiating A7 with respect to Yr yields 

0°76 Ll ay 1 ow 

ee 26S nee (AQ) 

ay? Gero lo ieee og 

Differentiation of A8 with respect to 6 yields 
Zz 2 
Be = eae) 
r9e* 


mee Climination Of derivations of ~ using AY and Al1O shows 


2 
¥ aa ee ae, (All) 
ra 0 4 


0° 6 
ae 








All is recognized as Laplace's equation in two dimensions 


ee) 


vo = 0 Gay) 
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From the principles of vector calculus it is known that 
in a conservative magnetic field region in space where 


mpm (, there exists a scalar potential such that 


me 76 (Al3) 


; Hy = a (Al4) 

Notice that the machine is modeled in two dimensions 
only and therefore H, = 0. Equations All or Al2 may now be 
solved with the appropriate boundary conditions and uSing 
Al4, the components of H may be found in cylindrical compon- 
Birt S . 

Using See of variables (see reference 1, pages 
434-435) equation All may be solved by first assuming a 


fetieion of the form 


ane) e= Rr) 0 (0) (Al5) 


berererentiating Al5 appropriately, 


(Al6) 


Q9] a> 
Dj 
i 
wW 
© 
a| a 
a" 
oO 
I 
wW 
© 
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Substituting Al6 into All 


ce OC += R' 4g += R 9" = 0 (Al6a) 


by 


And, upon separating variables 


(r*R" + ~R') = - = = k? (Al7) 


ee) te 


Where k? = separation constant 

Because the left hand member of Al7 is by hypothesis 
independent of pt athe the right hand member is independent of 
Meee tOllows that both sides must be independent of both 
r and 6 and equal to a separation constant, k*. With some 
miergnt as to the form of solution expected, designate k?’ 
as the constant. Equation Al7 implies the following set of 


Second order ordinary differential equations. 
meRe torR —- k°R= 0 (Al8) 
rk ©. = 10 (A19) 


The sign of k* is chosen at the stage of the analysis 
Bemrhat periodic functions which are continuous will appear 
in the solution of Al8 rather than exponentials (sinh) 
in the 0 direction. Equation Al18 is an equidimensional 


linear second order equation and may be solved as follows: 


oS 





assume R = a |e 0) (A20) 


Differentiating and substituting A20 into A18 


1) + ae a = = 0 


= tle (A21) 
and, 
a k —k 
R = ay. r+ Dy i | (A22) 
assume R = on k = 0 (A23) 


eeeerentiating A2Z3 and substituting into Ais’ 


ok tf rR! = 0 


peed 
jou 


i 
A ]e 


tim R= =—(inr) ' 


imeke 9= —lnr + in oo 


lea 


50 ote Bg bo inn; kK = 0 (A24) 
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Equation Al9 is an ordinary second order linear equation 


with constant coefficients, and yields the general solution: 


io. = C, COS k6 + dy Slice | Kk..2 0 (A25) 


O = Co =F d 8 , k = 0 (A26) 


Substituting A23 through A26 into A1l5 the general solution 


feel 1s found to be 


fen eae been (C+ D. In y)8 + 
2 eee + Br“) cos kO + (cr + By a) Sin k@] (A27) 


The following constant equalities are applicable to A277: 


A, =a 5 Ay = ay Cy 
De = DC By — Dy cy 
(A28) 
co 0°6 Ss oie 
Bor °o"o Py Sey 
In order that @(r,8) be single valued set as agen = qo = 


For the same reason, the periodic part of A27 must possess 
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a period of 27. This requirement leads to the permissible 


values of the separation constant k. 
eee 1 — L, 2, 3,422. Pep epee S. O.y Ste farce 


It will be later shown in this analysis that n may assume 

Only odd positive non-zero integer values. p is also later 
Meememed aS pole pairs in our analysis. If it is required that 
d(r,0) to be defined at y = 0, DS must be 0, and A27 becomes 

k ~k 


es) COS a um ah (Cyr + Dyer eisheee <p 


O(r,8) = T[(Ar< + B 


im 
k k 


(A29) 
Mees solution of A2Q subject to appropriate boundary 
conditions, must meet the following additional requirements: 
ine Be finite. 
ese a SOlution of All. 
Beecaticnry the circuital Or curl cOndition imposed 
by electromagnetic considerations. 
4. As will be shown later, @(r,8) must be an odd 
periodic function of 0 with a period of 27. 
Sensci@ering Only One particular harmonic solution of All, 


A29 may be rewritten as 


Weems 4 Be )eos kO+ (G-* + Dp “)sin ké (A30) 


me-alazing A, B, C, and D are functions of k. 
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Utilizing the technique developed in references 2 and 
3, A30 may be written in complex form, and the general 
expressions for the magnetic fields and inductances (see 
Appendix B) may be developed. The machine will be character- 
ized in terms of surface coefficients defined as functions 
of machine boundary conditions. These surface coefficients 
are then substituted into the general expressions for 
magnetic fields and inductances to find the specific solution 
desired. ThisS approach allows the designer more flexibility 
in the machine analysis than the more classic approach of 
fixing the boundary conditions and then solving the magnetic 
fmretd problem. 

Wiererag A330 in complex notation: 


k -3jk0 


® = RelAr™ A eee a 


Ea ] (A31) 


It may be shown that 


[to 
It 


Bete )D 
Allowing 6 to vary both in space and time, A3l1 becomes: 


Nee 7) Mb Calc) 


® = Re[(Ar® + Br‘), (A32) 
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fieeas now possible to examine the definitions of the surface 
@eetticients, S$, and the reflection coefficients, a. A more 
detailed derivation may be found in references 1 and 2. 


Using Al4 and A32 


ote lt Hi, (A33) 
when 

a Smuts ej ico? (Eo KY) y (A34) 

H = RTs) (Cae Pemeee res oa) (A35) 


A parameter referred to earlier as the surface coefficient 
which is the ratio of the radial to azimuthal component of 


mae t£teld 1s defined as 


s= (A36) 
6 
The complex surface coefficient is defined as 
H Ay me i ee 
a J | 7 (A37) 
Ho ayes is Be k-1 
Eoethe complex reflection coefficiant is defined as 
iit. S (R) 
i i eee (A38) 
_ 1 ~- 3 s(R) 
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and it may be shown that 








nae yes 
s(r) = jl —+ } (A39) 
1 + a (By 2k 


A393 allows s({r) to be transferred from one boundary to another 
Meee radial direction. The development of this expression 
1S addressed in more detail in reference 2. 

ieee bellow retlects the values of o« and s for the 


boundary conditions indicated. 


BOUNDARY CONDITIONS Ss QL 
Nothing Inside +7 0 
Nothing Outside — 7 00 
a O 0 +1 
i= > © co -] 


ition. SURFACE AND REFLECTION COEFFICIENTS 


Referring to Figure 1, field equations for the machine 
may now be developed assuming no magnetic material inside the 
driving shaft and no material outside the shield. At this 
point in the analysis no particular value of the shield 
permeability will be specified but will carry as the variable 
H,.- The air gap fields both inside and then outside of the 


field winding will now be developed followed by a development 
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of the fields within the field winding. The results are 
then extended to the armature and compensator windings. 
Based on the two dimensional schematic representation 
of the machine in Figure 1 a harmonic expression for the field 
current density may be developed. Figure Al depicts J (0) wise 


§ and is used as a basis for the following Fourier analysis. 











J, (8) 
imi 
SZ 1 
WY mw ows p 7 Owt 
p Z p 2. 
Ow Ze ite 2a 
a p 2 p 
Figure A-I| Field Winding Angular Dependence of Current Density 
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J. (9) may be represented in a complete Fourier series 


as developed below: 











i i __.2n6 _ 2n16 
J ¢ (9) Jo + (J, cos 5 + K, Sin ) (A40) 
n=] = = 
merce P 1S the period and equal to =. 
2m 
J = i ie F(G0)d0 = gi sP F(0)d0 = 0 (A41) 
fo 2) ai 
O O 
: 2m 
J - a 2 f F(®)cos eee deo =28 peas F(9)cos np0dd (A42) 
n p p TT 
=— Oo 2 O 
where F(9) = J, (8) 
Owf 1 Owf 2 
pI Z p 2 p 
1 i [ feos n p8d6 - fcos n p8de+ ff cos n pds) 
. 1 Owf 2n_Owf 
p 2 p 2 
(A43) 
Integrating A43 yields 
J _ Je [si AIDEN A ape 2 AAS ae gay ates Sane 
io i) 2 2 
-sin (2nn-2PS) | (A44) 


7 





Using a trigonometric identity it can be shown that 





AJ 
f .. npowf 
at sin 5 igh iereliel 
en i (A45) 
O n even 
We can show that 
ie 
a = ime ein oo” a6" = 0 (A46) 
Be PP» Pp 


If we allow Je to be a function of both time and space, A40 


becomes after substituting A4l, A44 and A46 


J COs (iit mous) (A47) 


a. 


li <4 78 


J-(8,t) = 
n 


in) 


where Jen are described by equation A45. 


The harmonic components of J_(9,t) may be written in a 


complex form as 


i, pe) = Relies Rites (A48) 


Where J = J fe Orme OCC ac defined in A45. 
ig) fn fn 


he 





Piewetiective field current density may be related to 


mie cerminal current as 





Mec: ft cif 
Je OR - 2 2 
E cee Are Gee ) 
Define 
a Rei 
6 
2N I 
J 5 re me a (A49) 
oa ay )Rey 
aid J becomes 
ela 
ON ¢ hes mp og 
; : sin —>y— mo oad 
nO -(1-y Reo 
Jen =- CAD) 
0 n even 


Referring to Figure 1 and references 1 and 2, the field 
components, H¢ and Hogs in each region of concern 
(Ro eis, °< Resi Res << Regi Reo S16 Ro;) may now be 
derived. The field equations will be treated as complex 
expressions, and the detailed derivation of these expressions 
may be founded in reference 2. 


Again referring to Figure 1 and reference 2, the field 


component in the radial direction inside the field winding is: 


We 








mguGeA = 2 Field Winding Pole Face 
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(A51) 


The separation constant initially introduced in equation 


Peers now defined as 


k = np (A52) 


Tt is assumed at this point in the analysis that the 
permeability of the shaft material approaches that of free 


Space, therefore a. now equals 0. This is done aS a matter 


@reproblem simplification realizing that Os could be carried 


throughout the analysis. A 5l becomes 


5 ey R 2t+k_p (2tk 
, dota 16 £O fe 
se oi ot) “Sa? | a) ne _ 
a 2 Gis its 5 2+ k 
EO: Si 
Zakk 2~-k 2-k 
eoe eG JES Bad 


— 


J (Sey 


12) 





In general, 


E. ) (wt-k@) 
H = Re LH ] (A54) 
therefore 
2+k 2+k 
Yin _k-1 Reo Rei 
Hep = Re | J- (+ —) | - 
nee 2 a R 4k 2+ék 
£O° Si 
Zhe 2-k 2-k 
Cfo “si (Reo ea ) acces 
et \& (A55) 
2-k 
Pmeorder to simplify A55, define 
2+k 2+k Die 2-k 2-k 
Sao) es Seat ay es ae 
QQo= ——— (A56)} 
2+ k 2-k 
and A55 becomes 
J 
Eno ok ; 
- ——— J (wt-k®0) AST 
Bet a te I~5 a 7 2k ) - (A957) 
ine mec 
pubstituting Jen = Jen + 10 Ianto Ad5S/ and reducing to the real 
Part only Be becomes 
y _ Je k-1 
rfi > (Ep) 2sin (wt-k6) (A58) 
a. RR.” 
iO Si 
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Substituting A50 into Ads, 


4Ni¢ T, eae . n pd 
re Qsin 


Heei = 7 ( 2k) y 


a 2 
Us Y oe CFOs 


Sin (wt-k0) 





(A59) 


A59 is an expression for only one harmonic component of Heese 


and the total field must be represented as 


7 00 AN, ¢ I, cok 
Heed a Rae 10 ie ae 5 a - ay) 
‘ te “ fo £O. Sa 
n odd 
KO ¢ 
sin 5 Sin (wt-k@) (A60) 





Relieving A60 of the requirement for time variation it 





becomes (W = 2nf = 0) 
2 oo AN, ¢ I. xl 
ae (l-=y2)R, 7 ear 
ae OBE Y fo fo si 
heoad 
k6 f 
sin = Sin k 0 (A61) 


Before proceeding with the field analysis an examination 


of a iemwmeorder. This derivation is patterned after that of 


reference 3 and utilizes Figure A-3. 


te 


Aeethe OulLer boundry, Beat S(Ro4 ) = -3, and at the 


Suter surface; 


ue 









Free Space 

#+ o 

Free Space 
HE o 


migure A- 3 schematic of Machine Shield Showing Parameters 
used to derive a,. 
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Wee so) eo {Reg ? 
U 
ee ui 
Beso i ue Ue ee ) 


+ - ie 
where aa and ae Bepreecehnue Incremental poOsltetons adjacent 


to R . It may therefore show that 


SO 
= LU 
7 Hey Ce ) 7 Hee Ce 
S(Ro, ) = = a) genie 
oleepe as ) 4- 
£0 * so_ Heg (Ro ) 
* 1S 
s(R.) = <2 s(R,.") 
u SO 
S 
on Ne 
and s(R )=ery— (AG6Z) 
—'"'so He 


Pemectomsrer S(R) across a radial increment use the 


equation A63 below '*? 


ck-1_gp2k-k-1 
ee a | 
_ * t LAR 


2k -k=1! ee 


Substituting A58 into A63 and simplifying yields 


1~ (By 2*—4s (R) (14 (2) 7%} 
s(@) = 30 | OOS (N64) 
1+ (=) pale Oks, [1~-(=) ] 


(ae, 





At the inner surface (r = eee, of the shield, the 


boundary conditions require 


SUR) (A65) 


Equation A6é4 is used to transfer the known boundary 








Seneition at r = aa (equation A62) across to aS Using 
equation A65 the desired boundary condition at r = rn may 
be found as 
Rene I R_. 
eer S1) 2k eee ay S1)¢K) 
= Ms Reo Ms R50 
s(R_. ) =-j - oo (A166) 
Oe are ey ae ie Rei, 2k 
1+ (———) + —— [{l1-( ane 
R U R 
| SO S SO 
Because S(Ro; )} = Seo (Roy) using equation A38 it can be 
shown that 
Ls eg Beh aa! 
fo ui J Se, Ro) 


A general expression for Oe, can be developed by substituting 


A66 into A67: 
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(A68) 





Note that eg Poenettyand in future «analysis the com 


plex notation will be dropped. Also note that for all 
windings a5 is equivalent and the first subscript may be 
dropped. 

In summary, the magnetic field density for the three 


particular regimes of interest may be written as: 


i RS eRe k = np not constrained; the most generalized 


ea! 


=eOrm: 


© AN, 1 ro k6 
He = bE —th ££ _  (E___) Qsin késin 


ied jatk nt (1-y*)9 
1 COOKS! 


wf 
2 





2 
WE £0 ee 
(A61) 
2+k 2+k as 


9 - fo —f1  _ 9 Si_ (R _ Rp. 27k) (A5S6) 
2 +k a ee 
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a defined by equation A68. 


























is) Rei k= 2, n= 1, p= 2 
4N i 
ee | 8) (0) Sin 2 6 Sin 6... (A659) 
rf (bows One RA. 9 edo Ra * we 
wf fo Osu 
Re, (= Ro." R 
ake a \ igo) 
= 7 OO Ros In(s >) (A70) 
fee 
Boe Hy oe 
ee ee a 
Hu SO S SO 
i = ee 
Ho Ned Ho Kes 
ae. et BG era) 
Pa SO S SO (A71) 
O 
Kot Ho Nef 4 
2 - (—-=)* + — [1 + (—==)*) 
u R U R 
1 - es SO S oe 
7 re (Sty 4 MO tg (SEY 
R u IR ! 
SO = SO 
allowing H_ to approach © reduces A68 to 
Cae ok 
O 
ier Resi He > ae not constrained 
00 4N 1 k-l k@ 
pe ee Se) 8 Sin k 6 sin— 
ae al) @ ni(l-y*)R. 7 R zh Z 
wt fo Si 
Nneeodd 
(AD) 


BZ 





Cate 24+ 2k 





(A73) 


6 


u 
on 





eS) Of Sim 20 (A74) 


fo 
= rt ee x eae 
Q n + Ro ln ( rm, CAW 5) 


Again, uSing results from reference 2 the azimuthal 
component of magnetic field density may be found as 


k-1 2k -k-1 
en ic +06. Ros ie R 2+k_» 2+k 


iO) fe 


era 


ee 
_eha 2. ae VAI R 2k 
Oo Sl fl al 


2k 22K eos 6 
R Ld ( R me e 
iy ee fo ia ; (A76) 
DAK 


again -.=0, and taking real part of A76 yields: 
5 pk-l 
— fn (——s—) 2 cos (wt-kd) (A77) 
H ‘ a aie 
Of1 2 Oo 





Oo. si 


Using A45 and realizing w = 0 it can be shown that 


5. ee k6 
Hines. — (a) en ee ee (A78) 
Oo 





2 


§ 3 





2 defined by equation A56 and OS by A68 


ie = R k = 2, n= 


fi’ 


ip 2s enor Constrained 


S 


a defined by equation A71l and 2 by A70 





25. . 
Hoes = ce ae WeOSte2ae Sin O of 
On cs 
fae R..; UP >: k not constrained 
eat S 
= -l; & defined by equation A56 
00 25 - ysl ke 
Hoey S = 23 a EeveOowih sl Sin 5 
n=l R_.2k 
Si 
mn oad 
i@ es Reji We + ©: k = 2, n=l, p= 2 
ee. -~1; 2 defined by equation A70 
~25 5 - 
H. . = (——) © cos 2 6 sin 8 ,. 
( TI R 4 WL 
“Si 








(779) 


(A80) 


(A81) 


The fields generated by currents in the field winding 


regi R Sg ck 5a see 
in the region £0 ae 


reference 2 


will now be examined. 
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From 





H = j — (ee i [ fo sD ae 2k 
inne, 2 2k OK eee Fam 
Ua Gee Roe 2+ik ——- 
a 
] (A82) 
2-k 


jen Ue 


Poemng the real part of A54 and considering all harmonics 


Soouning nothing inside the driving shaft and a 





yields 
pe j Pe Sas R (2k =k R 2+k _p 2tk 
i 8 O Sl fo i: 

Hepg = t 8 Se [Se] (2) 

n=1 Rime Ee 2+k 

oe Si 
meodc (183) 
Sin vk 6 Sain ae 
2 

On defined by equation A68. 
Reg ue Roji k=2, n=l, p=2; He not constrained 


Qe defined by equation A71 
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R Weer ier kK Not Constrained: u > © 
© oli ha Ue 





00 NP reeaoeege TKR otk 
Heeg = t 2 oe (eR) (SS 
n=l Res a 
Moda 
KO ie 
Sin k Osin 5 (A85) 





i 
rfo nit oe) so sin 20 sin, 


(A86) 


The azimuthal component of the field outside the field 


winding is now found from reference 2: 


J k-L 2k -k-l 2+k 24+k 
cn hl PO ts ie | £0 Ry. 
2 pte 2k 
O.OR. -Q._. ; 
oO: Sy: fl al 


H = 





mae Se? 
Cy (87) 
ane 
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des = 0, taking the real part of A54, and considering 


all harmonics yields: 














bs ki oe R 2k_-k-1 R 2+k_p 2+k 
ae fi Oo Sl EO, ah 
cee Unt CCC ) 
n=1 a oRo: 2+ék 
nOodd 
k@ F 
Sain —;—cos k @ (A88) 
Ol defined by equation A68 
Re ae ec te 2, te not constrained 
O si S 
PR Ce Ror! Roe Ro 
H a £ (ee Sy te fal VCOs 9-2 -oSin.w 
8fo TT 4 4 wt 
Oak © 
O Sl 
(A89) 
Oe defined by equation A71l 
R Seek .- fk not Constrained, uU. -~- © 
fo sl Ss 
a =-l 
O 
a a5 yk-l R (2k,-k-1 R 2+k_p 2+k 
1 Sl £© jae 
fot = 7 a an S 2k ) ( es 
pe si 
no.oad 
Ke 
Gos) ko sin 5 (A90) 
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Reg eis: ~< Roi = aie a= 2 - He -> 00 





a” ae a 
[oe GOcmo. 0 Sin. © 
wt 


(A91) 


The field expressions within the field winding are now 


developed. (R 


ci ie ie be Reo). From reference 2. 
Jen pK-1_y R (2k -k-1 rkt2_p (2+k 
7 4 = (_— On ee \f dea ae (2k 
—rf Z ae (2k y R (2k 2+k fae 
OS 1: ig isla 
k-1 2k _~k-l 
(r? Kor, Ky Jen ie “Oe Ros -r 
= er ae ) 
2 k Z - oan R (2k 
Si iL. ad 
Reo S 2+k oe a Pe a 
[ Slee (Aa2) 
26k 2-k 
des = 0, taking real part of A54 and considering all 


harmonics yields: 





e oJ pk-l_y R 2k ~k-] 2t+k_» 2+k 
i O Sl : lem 
Ho = » n ( Ik ) ( ) 
n=1 vas aN 
n odd 
OAs 2-k 2-k 
F ie to sili 2+k - O16 ae (Reo “Rei f 
a R ans 2+ ék 2-k 
O Sl 
k@ ¢ 
sin = sin k 0 (A93) 
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OO defined by equation A68 


Ry; << Kei R=2, n=l, p=2; He not constrained 


2 ig OO Ney Te a : 
eee yf ) ( ) + ( ) 
OL : 4 Oe Ro? 
SUL oy eal 
Reg "FE" ae 
= + iv : : oO 
( 7 le mm = es ales! 0 poin 270 (A94) 


Qo defined by equation A71 


Loewe k= - | + &- k not constrained 
fl fo S 








zs og k-lp 2k -k-1 2t+k 2 (2+k 
ae = | een ee eee 
ie fe al ni R 2k > tk 
S71 
n odd 
K-11 R 2+k _2tk R (2K ip 2-K ek, 
ho Sl LHe: 
+) | + 
2+ ik mn. 
Sl 
KO ¢ 
Sin 5 Solace! (A95) 
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ES < ° = = pe. 00 
Rey r Reo? k=2, n=l, p=2; We m2 








LT u 3 yo 4 GH 4 
_ ale aes ie - Sek r ae a 
a (+ (——) 
Re 4 R.; 4 
4 es 
+ —— 
Ro; Ti = )] sin ae Sin ae (A96) 
Again, from reference 2: 
k-1 2k -k-1 
Jen ie ta, Ro. ret k_p 2tk 
we | 3k ye = 
—_ CG. Ru. -Qa,_.R.. : 2+ ik 
O Sil fl al 
2k, 2-k . 2-k oy ee ae 
ere Ge ~R.. ) J fy an 
ee sie 4 fn; — 
ee Nek AE a en 
Oo sl fi al 
Resi cr ” (Re Pee) 
[—-—_-__—-. ——. - a Ros ] (A97) 
2+k 2-k 
de, = 0, taking real part of A54, and considering all 
harmonics yields: 
oe k-1 ro eee ek 
Hos — (~—s,) (=) + ( Dk? 
=) - aR. Ya ye a R 
Oo Sil Oo si 
n odd 
R 2+k_ .2+k R (2K ip 2-k_2-k, 
te it Si ie 
etal p atk 
k@ F 
sin —;— cos ve (A98) 


a0 





Oo defined by equation A68. 


3S < . = = ae i 
Ry; iis Regi R=2, n=l, p=2; He not constrained 





2 ee See | 
aR ele = Jasin oe cos 2 90 (A99) 


Sah 5 


OS defined by equation A7l. 


iru < : 2 
Res Ie Reg R not constrained, We 5 








co 2J ois eee J ek 
H ae SD —--- (——_.-) (—_—__—__) + ( ) 
a n=.” | 8 2+k ao 

Sl Sl 

n odd 
R 2t+k _.2tk R 2K ip 2-k_ 27k) 
fo Sl fo 
a ee ee eee 
2+ k 2- k 
k@ F 
Sin = cos k 6 (A100) 


oa 





cee = k » k=2, n=l, p=2? Uw. > @ 








Pal rt-=-R oe R Hoyt R 
1F il - 
Hog = oq L() (4) + () (42 — + Sn £2)] 
a 4 R 4 
Sak Sl 
sin Of Sioa @ (A101) 


Referring to Figure 1, the field equations as developed 
above are valid for the fields generated by the current in the 
compensator if we replace 9 with 0-0. a is defined as the 
angle between the direct axis field component of the field 
Winding and the direct axis field component of the compensator 
Winding. We may also find the field expressions for the 
armature by making the appropriate substitutions as outlined 


in Table A3. 


a2 





Substitute term in column B for column C to find fields due 


meecurrents in column A. 


A . B C 
Compensator 0 
WC oe 
a J 
Me ALS 
Be te 
A arn 
eo EO 
Z y 
6-a 0 
Phase A or Armature O+¢ © 
ae wf 
es J 
Mee ee 
He Ie 
Sa SH 
ae Nae 
x y 
27 
Phase B of Armature 6+o+ ae s) 
4 1 
Phase C of Armature O6+o+ > 6 


TABLE A2. Substitutions for Field Expressions 
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SEPENDIX 2 


MACHINE INDUCTANCES 


iiieerlux-current relationships for the machine may be 


expressed in a matrix format as 


a Pwo ioMee le 86. Te a 
ro Oe ee: oun OE Eo, 
oie 13 3 864s Boe Mee oe Fe a 3 
MA aie "Os ee Le er ae 
OG et OF a oe fag Ue | Me 


In the above expression, for clarification a, b, and c 
(armature phase subscripts) have been replaced with 1, 2, and 
Beerespectfully. 

In general, the self-inductance of the field may be 


found as 


heterring to Figure Bl the £lux linkage may be interpxreted 


memcnat £lux, © passing through the coil. More exactly 


X = N,9o (CBAs) 
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Saiesingie Turn Coil Flux Linkage 


Figure 





where for a single turn coil N= and 


me / B+ n da (B2) 


equation Bl becomes 


me! Ben da (B3) 
S 
By referring to Figure B-2 this concept may be extended 
to determine the flux linked by a thick winding composed of 
many turns. The differential number of turns in the 


differential element rdwdr may be expressed as 


2 he 
2N, 
EN, 3 ps ee eg 7) (B5) 
Gee. °-R...*) 
W wo WL 
where dA = rdrdyw 


The flux linked by this element is given by 


coe od NL | (B6) 
2N } 
d*} = ————-___—___ dA (B7) 
Sei OR) 
Ww wo wi 
where ® is expressed in equation B2. If we are determining 


the self-inductance of the winding, ¢ is generated by currents 


mlowing in the winding in question. Otherwise (mutual) ® is 


a 











Memerated by some other current within the machine. Within 


the air-gap we assume the fields are conservative and may 


feeerecsSS B as: 
B=u_H (B8) 
and, B/ becomes 
d7A = en MoH - n daj a’N, = (re u_H «= n da] 
(B7) 


Treating the field winding self inductance, equation B7 


may be written as 





wf R 
P 2 ae v 2N, -*rdrdy 
L-== f/f i bo eH (ee, oO) ead Z | 
if I. 9 Ort @ _(R..2=-R..2) 
_ wt R jeu wf fo fa 
2 ies! p 
(B7) 


Using equation A93 we may integrate B7 and obtain the field 


self inductance: 


eoe 
fore) ee a eee 
besin® (>) Nie * slo 2-k ee cee 
Le = 1s Brn eee) Os 2 lane a ) ae =) 
n=l aon pie. (k“-4) (1-y*) S1 
n odd 
to [(k-2) + aa y(24K)1 | (B9) 


oo 





when OS is given by equation A68. 


Let We oO, k not constrained 





a) — | 
© 
ma er wf 2 
i abillde ai Say ae alae age (2 (KZ) (qey2tky 2 
© n=l n®pme, 2 (k2-4) (ley2)? Kt? 
MwOaad Mi 
ee Oe 2+k 
(am) + (ke 2) + ay = y* (24K) ] (B10) 
Sr 


This result agrees with the tabulated value for Ly in 


reference 6, and can be obtained directly from equation 60 of 


reference 2. 


k=2, n=l, p=2, He not constrained 








re 2 
3 aoe WE) NEE Seg Dil-y ye ($f) 4 
it 2 ee 8 Roce 
OT WE ae = 
R 
4 l-y* | (l-y")? Zone 
ey ny a oe ] (B11) 


On defined by equation A7l 


9 





O 
ae 3 2 
Meee et ENS Ci-y")? Feo, , tay" 
be = te ee (F8)  ytan y + 
10 ¢ lice yaar) Sl 
(B12) 


Using the same technique discussed above, we may find 
all machine self and mutual inductances, allowing for the 
angles between ranous windings when developing the expressions 
for mutual inductances. The results of this analysis are 


listed in Table B-l. 
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die ev, S) SAB tod E 
MACHINE INDUCTANCES 


See Appendix B for expressions for L;- 











k = np not constrained, H NOt COnNStEralned 
2 . Wa 2 
7 a eee 2 Nes ae 2-k Dee G 

L. = 2 —— 2 (Pym) (1-x ) 

n=l a.n®*p 0 *(k*-4) (1-x?)? 

Me oad 

Rao, 2k 2+k 
Ga Rae “=x (24k) ] (B13) 
Sal | 
OS given by equation A68 
k = np not constrained, He eae 
: oo ; 
“ i i meee Z [ee * Mo es Ze A@ wa 

BoB EE 2 5) (1)? CE) 

noe nN pT Fag) Cl es) ai 

ig fexellel w 

+ (k-2) + 4x7t® ~ x (24k) ] (B14) 

Co) aamaee 
O 
7, 0, p=2Z, Lt. > © 
ce = —) 
O 

Somme GG N21 De IR oa 
L = —___Wa 3 2 0 (Lire) (282) "4xhan eas = ) (B15) 

Tey le) Si 


JOH 





TABLE B-l continued 


[ae n-t, P—-2, U_ not constrained 


ee 


le 2 
0 ; 
ee wal "ta “ato = [_(1-x")? Pao,, 
a Drage oO i 
see (l-x*) Si 
F Ro 
4 (Ciara ene So) 4 
+a [x In x + n eras eam (==) ] (BEG) 
Roi 
Qo given by equation A71 
k = np not constrained, H not constrained 
“we 
oO Aa 
a l6sin® oe De Tells 2-k 2+k 
L. = 25) tt 2 (yap) (1-2 ) 
Dele aon eg 
ig Tejskel 
Fo, 2] 2 
(==) — a, [(k-2) +42 TS 2” (24k) ] (B17) 
Si 
Og given by equation A68 
ife—otip Not constrained, Woe 
eae : ee : 
. a) eel 2+k, 2 
L, = 2 a [2 (S54) (1- Z ) 
m= nape. (k*=4) (1-z*)* 
WC 
ey erate! 
Nee 2+k 4 
=a + (k-2) + 42 + 2° (2+k)] (B18) 
Si 


LZ 











PABES B-l1 continued 
f=2, n=l, p=2, He NOt CONnStrainea 
2 2 
0 y 4 
= Bsa we) Ste *eHo | TG alone 
c aa eee ee & eee 
Sa Ue {1L=Z~) S1 
t R 
, (laze. Gaz) 1. Cov, 
eee eg eg ee) | (B19) 
Sa: 
a given by equation A7l. 
k=2, n=l, p=2, Ue ane. 
meee 2 
2 Sees Se ee wee (l-z*)? ic \ ty 1-z" 
La >), ar mmei aa ap VA Mol yA as zi ] 
eee (1-2°) Si 
(B20) 
a =-l 
O 
k-ap not constrained, He not constrained 
ale, 
fs a, wa > 2nt 
l6ésin* ( 5 Ny, hus cos (~,—) ok Dee 
L b » 2 (rap) (1l-x ) 
eeoem=t o n®* po _°(k*=4) (1-=x*)* 
O wa 
neodd 
i + a [(k-2) 44x97 Sax (24K) J (B21) 
Si 


Qo given by equation A68. 
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TABLE B-l continued 


k=np not constrained, He ee 


k0 
a a2 wa 2 ¢ 2nT 
L ee 5 pas ( 2 ee ae IE 3 a) ann? hye 
ee y= 1 Be i (ena) (lima) so 
page re Ke | 
6 2+k 
Cas + (k-2) + 4x - x (2+k)] (B22) 
Si 
a= —1 
O 
k=2, n=l, p=2, He not constrained 
eae 2 Zu 
L a 8sin Se ae ue SSIS) (l-x") AO) 4 
eae Gg 56 2 (1-x2)2 Ped 
O wa 
18 
4 “ese 7 (ules) AO) 4 
+a [x*in x + =e 4 Se (ae) 4) (B23) 
Sl 
Oo given by equation A7l. 
k=2, n=l, p=2, ais 
Oe = = 
- 2 2 Zit 
: _ 8sin Cee Ne Ru, cos (=) (ex)? (£80) s 
ab Tere (la) 8 ee 
ae Di: 
+ x*lin x + a ] (B24) 
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TABLE B-l continued 


k=np not constrained, He not constrained 





2 ae 2nT 
2 00 he (oe = asec 
7 16 ae Nea Sale 5 ) (l-cos 3 ) 
L “Lip So y CC 
a ao TO eae jie n=1 n? snp 
Oo wa 
neoad 
(B25) 
mee 2 (2%) (1 - 24K) 2 (80) 244 fC ee = On 
snp Kee e ie Or : - 
(B26) 
O defined by equation A68. 
k=np not constrained, He + © 
k8 
eo uN. 00 sin? ( 24) (licesseey 
= ao ta 2 3 
L -L oo ns a con 
— * Toei Da) n=1 i 2 
n odd 
(B27) 
k-2 Pare Rao 2 2tk 4 
— ——— _ =) —-2)+ —- +k 
en 2 (Fy5) (1 x ) ‘Rage (k-2)+4x me eels) 
(B28) 
Q = -l 


OS 





TABLE B-l continued 


k=1, p=2, n=l, D not constrained 


ole Ne 





aC. ta eee Peal) 
iL SIF = ee Sin’ (9 _) (1 -cos—) C 
a ab at6 2 (1-x?)? ie : ar 
Oo wa 
R 
7 _(1-x")? AO. 4h \ (l-x") 
Co a 5 same + a x lige 2 ae aa ee + 
Si 
Oh Goven by equation A7l. 
il, D=2, n=l, H. 7 oe 
Qo = -) 
O 
Seen ioee Nee 
bo bap = #2 =e — sin? (6,4) (1-cosS) “52 
10  *(1-x*)* 
wa 
bs (a 4 BO ar b 1-x" 
Coo = 7 ae ee <toleuipores ae r 


PmMpenet COnStrained, We not constrained 





Moe = 2 Mot Cos To 


n odd 
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(B20) 


Ob 0), 


(Bal) 


(B32) 


(B33) 


/ 


TABLE B-l continued 


where 


s 


kO k0 














-} Zs 
eo tay ") Lm, NH sin(-—"")sin(—™*) R 
ae = - ie 2 2 “£0, k 
mi (RE CnpEt 
am ae Heer Le Vee Ss: ao 7 
(B34) 
R 
ee = ae +k AO ie 
F 7 [di a (Soe) (1- x? ) (ep) 
£ Se 
ripe ; (335) 
k(4 - k*) 
a given by equation A68. 
k=np not constrained, Me = 
Mog = a Mf cos k @, Qe = -] (1855) 3) 
n=l] 
n odd 
k@ k@ 
pesca Won: wf SoHE k 
32(l-y ) imu N,N + 5in ( 5 Saba 5 ) (5—) Sort 
M = ad 
ge 
2 
mn “8 (l-y*) (1-x?) 
wa oe (B36) 
ar 2-k 2tk, Rao. 2k 
(ul ar (s5;,) (- Xx Se } 
C 7 eS SS (B37) 
npt k(4 - k?) 
k=2, n=l, p=2, Me not constrained 
: —-321mu alles Nig Sin(0 4) Sin(8, .) Re ; 
M ay x UES en Sa a LL (2) (l-y‘)c 
Z2E 2f 
a 10 z 6 og tl- y*) (1l-x*) Rao 
2.” (B38) 
R 
oes - (1-x") (_ao)s 
Cor == 3 In x a 35 (= ) (B39) 


Si 


dO 





TABLE B-1l continued 


OO given by equation A7l 


k=2, n=l, p=2, i > © 


¥ 321m ON, Nig sin(0, _)sin(® .) Reg , ' 
M DMM =-—-—-ssa. > (-——) (i-y') C 
= 70 6 ,(1-y*) (1-x?) Rao 
wa wf 
R 
7 Cea eacnt 
wap og 1 * + 3D ee 


k=np not constrained, He not constrained 





M = 2 M cos k (¢6+a) 
ac ne 
n=l 
igi TONGUE! 
where, Ko poe 
: : ~321mN, Ny usin ( 5 ) sin (—s-—) | coy k 
1S O m6 g (eed: ) as npc 
O wa wc 
R 
2-k 2-k 2+k ao, 2k 
eos Re 
C a Sl 
mpc Nef alge 


6 given by equation A68. 


k=np not constrained, Me oe aS 





108 


(B40) 


(B41) 


(B42) 


(B43) 


(B44) 


| 
' 





TABLE B-l continued 





Beh ad caer 
Se i a ee 
” 7 321mMy ON, NL Sin ( 5 eae } Bco) k 
ike. 2 Be ae IR NPC 
™7 Oe On (l-x*) (1-z*) aO 
i ot SO 24k, (a0) 2k 
pug ey 
C = Si 
nes fe nes 
i-2, n=l, p=2, M not constrained 
7 : ~321mu ON, WN, sin(® .)sin(9 oo) , CO) 2 . 
Ae peas ee R ae 
ee ree ae, CL 2 2) ao 
R 
2 7 (1-x") aod, 4 
eee fe Oo, 3 CR 
Si 
ab given by equation A7l1. 
k=2, n=l, p=2, He + 
me 1 
O 
, 7 321my N,N, Sint? )sints Reo) 2 z 
a oe (x) ( 1-2") Rao Bae 
wa WC 
R 
eel: (1-x*) ao, 4 
Con ars oe 6 Se 30 ee 


IEG). 


(B45) 


(B46) 


(B47) 


(B48) 


(B49) 


(B50) 





TABLE B-l continued 


k=np not constrained, He not constrained 











Men = M fe cos k a (Bol) 
= 
meodd 
where, k 86 k0 
2+k wf ; wc 
P 7 Bale yy, ) ImN, Ny eH sin (—5—) sin ( 5 ! co 2k. 
nfc 2 Make 5 este R HpLe 
ain 8 ic we it y~) (1-z*) fo 
(B52) 
Fyne? ® = og ins Pe Wee LS 
ey ods ee 
iO = B53 
Bec rae kk) ait 
Oo given by equation A68. 
k=np not constrained, Wey lea 
M = }, M cos ka, a = -l (B54) 
PG: _ ec O 
n=l 
n odd 
k0 k8 
2+k : wf WC 
Sly )ImN, Niel, sin (—5—) sin ( 5 ) Rook 
snéc ~ m20 0 _(l-y?) (1-z?) ‘he “npfc 
- we wf Y 
(B55) 
Daa 2 =k De) Nae, De 
= Z, + (sa}) (a2 \is—) 
C _ Sl 
iit) Die: eae (B56) 
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TABLE B-L1 


Conernued 


n=l, p=2, He HOw GCOnstrarnea 








: ~321mu (l-y )ImN, peoind sin’ Roo ; 
ar Gee) lees Re “2g 
OF ve wl x : 
(Bo 7) 
IRE 
ess 7 (1-z") CO, 4 
eres 8 PZ Gap SR? aD) 
Si 
OS given by equation A7l1. 
=2, n=l, p=2, ie > © 
r 7 321mu ON, ON, -sind Sin® mec Ae — 
Dae R Y Zine 
ao 
(B59) 
R 
(1-z ) CO, 4 
= ot B60 
CoE 3 In zt 35 Re ( ) 


Lda 





We may determine all additional machine self and mutual 
inductances by making the substitutions listed below in 


Table B-2. 


Te Obtain Use Formula Fgquation Nos. Replace With 


ie ag 
M M B33 through $ (mes ou 
ong af B4l 3 
M M Bere clolig@sbisile w) db + ae 
Cl af B4l 3 
M M B42 through % ee 
be ac we ? eae 3 
B50 
21 
M M B42 through d+ a ¢ +0a+ — 
Ec ac B50 5 


TABLE B-2. Substitutions for Machine Inductances 


eZ 





APPENDIX C 


MACHINE INTERNAL TORQUE EXPRESSIONS 


Using the basic concepts developed in reference 7 and 
Other sources, excressions for machine internal torque are 
developed. Specifically, expressions for torque on the 
armature and for null torque on the field with a compensating 
feng Glther within or outside of the field are presented. 

In general, the Lorentz force is representative of the 
total magnetic force on a charge of moving with a velocity 


V and may be expressed as: 


fe oV x B (er) 


The force density, F, in newtons per cubic meter can be 


mecainea from Cl as: 





ees eee (C2) 


al 
I 
_ 
+ 
S 
I 


6 v>+0 Cae) at 
F.q., and V. refer to all particles in é6v and B, is the flux 
a 1 ald . ae 
density experienced by qa: Assuming all particles in $v 
experience the same flux density B, the definition of free 


Burrent density and equation C2 yields: 


F=0d0.xB (eS) 


Lae 





Equation C3 may be uSed to express the average of forces 
Simetne Charges flowing in the winding. This is equivalent to 
the force on the winding 1f there is some mechanism by which 
each charge transmits the Lorentz force to the material. 
Assuming that the charges in the winding conductors behave as 
particles moving through a viscous material, the force acting 
on each charge 1s transmitted to the medium by the viscous 
retarding force, and equation C3 represents the force density 
experienced by the medium (7). 

It is also assumed that the windings have a constant 


permeability approximating that of free space and that 
f= HoH is valid within the winding volume. (C4) 
Substituting for B, equation C3 becomes 
x. -H) els), 
Expressions for radial and azimuthal components of H are 
developed for various boundary conditions at the shield inner 


radius in Appendix A. 


Equation C5 may be expressed in a matrix form as: 


1. 1, tS 
ip = Mo 0 0 ‘i; (C6) 
i H, 0 
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or 


F = Hod (HL. Eee. ol.) (C7) 


Meee Wanding Null Toroue Expression: 


memeerring to Pigure C-1, the torque on the field winding 


may be expressed as 


7 = wee <P) dy = Jes x F) rdrdddz (C8) 
p ae 
= (ees Ed EC. 6 (C9) 
i Le 
ii. ly 1. 
ao Fr = r 0 0 = <F, 1 (C10) 
Z, 
Eas Fy 0 
Using €10, equation C9 becomes: 
T 6 _ 
ae f f r-Fedra0. | (Cig) 
g Zz 
2 Be RAs 
ir 3 


Equation Cll represents torque only on one pole face, and for 


mee entire machine, 


A (5° Z a 1) 
eee, ie SD ® | 
i ay Me Jf Hod HY Gras 1. ( 

i 9 Res 


> 





N14 9WOS | 


99D4 a/0g Burpuim play 


|-9 andi 3 





eG 





Consider the machine configured as in Figure 1 with the 
compensator located between the armature and the field. Using 


Bauation A&83, ie may be expressed as 





= 3 : 
HL a. 5 cos k (8 + o) Be + hee CG le3) 
or, 
a 2 agg + si 
‘f : 2 = tea oR k k 7 3g Rao. KR? : 
ig aan ni A. aR 2k 2 2+k 
; On Sil 
Meoad 
ae Rog - 
Sin (—5—) cos k(@0 + ¢) + JA atk ) 
k6 
sin (=) Sin k (60)] (C14) 
betsane 
aS Rao "Rai Soe 
a ae ) sin(——) cos k (8 - 9) 
2 
Rog “Rey eer: 
4- J ( ne ——-) Sin (—;——) Sin k (6 - a) Ces) 


Whereas, equation Cl2 becomes 


00 4pu_J cou oreo 
Be OS ee) ade 
f n=] ae: | a R_.2k 
n odd (elG) 


[ei 


~ 
Bo 


WL 





Mmeeagoating Cl6 and dropping the vector notation, C16 


becomes 
_ Zien 2+k | 2k Aa a 
i | uel: oe Reg AG RG (Reo ie 
A a nv 24k 2~- k 
im Odd 
ii 
Sl 
oe OC, tk, 24k Pee. ee, 2 
nel Fo ee! ee Zk IO neat O a 
Zz 2¢k Csi Dom Kk 
nao Kk 
Cus 


and integrating equation Cl6 across the pole face 








i aia 
eo > n ix 
he mS ; [ C,COS k(0+$)+¢,s1in k (8-a)] dé 
nodd — wf 
2 
where 
ous = — He 
ese 
ae ie kO 
eee 2k 
Mee integration of C19 yields: 
00 n 0 S) 2 
oe =) wel sin. i (e+) -sin k (o-—=) ] 
g k 1 2 2 
ie n=l] 
ie Oad 
a 0 
“C5 [cos ie Wea) -cos k (Sta) 1 


ME Es: 


i 


OC L7) 


(Cis) 


(Cio) 


(E20) 


Kez) 


(on) 





Mach reduces to: 


oo kO k@ 
an : Z | 
=e call ee eosin Cease ai on 6 sin (—-“4) cos Ko Ce 2S) 
# k 1 eZ 2 2 
L n=l 
n odd 





Pee quation C23 the current in the armature winding is 





Specified to vary as cos wt. Equation C23 may be expressed 
as 
— = o.cos kd ~ o.cos ka (C24) 
i a 1 y 
n=l 
ne Odd 
when 
k0 
2n wt 
me & ost 
1 
2n re : 
eee 2 


The condition for null torque on the field winding becomes 


Org) es: 


Oo.cos kd - o,cos ka = 0 (C26) 
a # : 


mo ac 
Which leads to a set of null torque requirements encompasing 
Pinodd harmonics. 

Considering only oe first Hermoenive, equation C26 yields 
as the condition for null field torque. 


0, cos ko = 0.,CcOs koe (E27) 


ode 


ana cos k $¢ = J A ,cos k a (eZ5} 


JO IGS. 








a. v 6 
where D4 .4k . k ae 
is -R, Sin ( ) 
i 2 ( CO ent ) Z C29 
me 35 I4k Sew KO (C29) 
a R -R , é wa 
ao al Sin (—y—) 


For a four pole machine, considering the first harmonic 


Only, equations C28 and C29 may be expressed as: 


Ses 2 0 = Joh elekey eel (Ca) 
Ro =k, sind 
fe = 2 co al WC 
i 3d, R so Fr sin wa 
ao ai 


Consider the machine configured as in Figure 1] with the 


compensator located between the main field winding and the 








earero. USing equation A83, He may be expressed as: 
H =2cos k(0+ 6) H. +H... (C31) 
r 2 rao ee 1 
- f 35, phot Rei ~k-1 as i 
He = 2 am lee ; ) ( 7k ) 
n=l] QL 
, si 
il @reUG 
eo coo Koo 
a icos k(e +o) 4a) @ sin (-—"—) sin k (0-a)) 
2 Coa R Zine ce Z 
Oo Si 


(C323 
where 2 is expressed in equation A3, modified by Table A3. 
PWeSserevting C32 into Cll and integrating yields: 


cos k ¢ = we A. cos k a (E33) 


110. 





where 


sae 24k 24K 2+k +k 
= Ys Sin ( 2 (Reo ae: ) (2-k) (R O on ) 
me 3 KO saa cat ve err 74 
ae wa ey +k at = oe ae 
Sin ( 5 ) aie Ro ) i al 
-a (2+k)R aor age ae 
O Sl CO alk 
0. ee eS 2, | (C34) 
O i ge a ine feo 


FOr a four pole machine, considering only the first harmonic C33 


and C34 become: 





cos2¢ = J, A, Siem dees) 
2 erie (Ree a pee Roe AG ke 
ie. = WC in® a Ge Ga: O SL 
2 3 b \ ee ie \ 
1 ee cee ee ail faa rates aes ope ae 
ee 
cae 
eke) (C36) 
jee 
ase) 
ea 


Note that for most windings where z and y are greater than 


a A, meances for the four pole fLilrst harmonic case to 


approximately the same value as Aj, if we consider a laminated 


aron shell. 


Armature Torque Expression 


: 23S] mee reaton ih « “— © <= R is 
The expression for Ht in the g an BS 


teal 





frome tt -. + Hs (Cae) 


USing equations A6l and A61, as modified by Table A3, 


equation C36 may be expressed as 








00 5 ee kO 
» ; | | 
HS iz, oa - ay) [5-8 _-Sin 5 Sin k 0 
odd Oe Sa 
ko 
pc ecin = Sin k(0-a)] (C37) 


where aoe te and 2 are expressed by equations A68, A56, and 
A56 modified by Table A3 respectfully. 


Peeng equation Cil and integrating through the armature 


winding, the armature torque per unit length is 


24k 2tk 











Ai ie Vian R Ro. k 6 
a Oa | ao al ; wa 
ae 8 Sa > 
a n*T Oh eas 
emeul 
KO of k6 a 
[J ,-Q-sin 5—Sin Say ss J a sin 5 Sin k (o-a) ] (C38) 


Mien K=2, n=l, p=2; equation 38 becomes 


T. 6u J sind. ; ; 
— u en : 8 s 
Tan manne eG Rai ) {J -&-Sind sin 2 
a IG, tl. 
© Gu: 
+ J 2 _ siné Sine 2 oo) (e232) 
Gc wc 


eZ 









When 2, and 2 are Crprecsed by equations A56 and A56 modified 


im 
by Table A3 respectfully. 
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APPENDIX D 


MACHINE SiZiING AND LOSSES 


As discussed in preceeding sections, this thesis 
addresses machines which may find applications in SWATH, 
iyarOrLOl1l, and SES vehicles. Current designs indicate 
machine dimensions and power requirements as indicated in 
Table D-l. Note the relatively large 2/D required for 
Meorotoi)] applications. 

Using the shaft horsepower requirements for various 
applications as listed in Table D-l, an estimate of minimum 
shaft radius may be made. This is equivalent to estimating 


the machine armature inner radius, Roi 


Shaft shear stress Tg, May be expressed as (see reference 
8 ) 
UG Rae 
Ea = a (ole 
i. Zz 
PWR TR 
where My = 7a ; I. =. Seen! assuming a Solid machine shaft. 
m 
(Dzy 
D3) 


P@estrituting D2 and D3 into Dl and solving for Rog: 


PWR(60))1/3 (D3) 


g 
TI To UN 


aa a 


Where N = machine RPM. Table D-2 is a summary of machine 


Mamimum R_. for aT Gm sO pO OUN pe Ik 
al 02 
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NAVAL NUMBER Sor” APPROX. PE PROX 2 EKO, 


APPLICATION MACHINES MACHINI: MACHINE MACHINE 
REOULRED ee LENGTH SUAFT IIORSI= 
ALLOWED ALLOWED & POWER RATING 
(in) (yy) D (each) 
D 
SWATH 2 96 120 ee 40,000 
Woee, Fig.2) 
iy DROFOILL 2 Sil FOO tea 20,000 
(see Fig. 4) 
Sis 2 V2 90 | ees, 30,000 


(see Fig. 3) 


TABLE D-1 Current Machine Vehicle Application Parameters 


Minimuro 
PWR N Roi Possible 
(ho) (RPM) Cant Sopp lication 
24,9) ONO) 1800 LES HYDROFOIL 
30.000 900 a6 SES 
nO OU 180 Grant SWATH 


TABLE D-2 Minimum Armature Shaft Inner Radaii for Indicated 
Naval Applications 


25 





SS 


I OS — 
= 


mre ttachine RPM is a function of the drive train 
Beeeecteristics between the motor and the propeller and 
propeller characteristics in general. Machine RPM may be 
adjusted by varying either the number of machine poles or 
Mmeequency Of armature current. The shaft radius actually 
S@eweyed 15 increased to account for some margin of safety. 

fn examination of Figure 11, reference 7 (approximately 
meeeecteed herein as Figure D-1), indicates the approximate 
domain Of magnetic shear stress to be expected for machines 
Beyer tOus power ratings. As a design guide the machine in 
Sieemeresentation should fall within the cross-hatched area 
indicated in Figure D-l. 

The magnetic shear stress may be expressed as 


Pe (PWR)P (D4) 


2 
(27K, 4) Qf 


Figure D-1l and equation D4 together may be used to select a 
"reasonable" value of Onn for any particular design. Equation 
D3, D4 and machine power as expressed by equation 12 are 
used to formulate the nachine physical size exclusive of the 
shield. 

Assuming an iron snield is employed, its thickness may 
be found by integrating the magnetic field density across the 
pole face of each winding at Rowe Summing the resulting values, 
and determining the thickness from an assumed maximum magnetic 
flux density in the shield. AS a Gesign constraint, a value 


pieeoeteslas 1S assumed to be the maximum value. 
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Referring to Figure D2 and assuming that all flux from 


Gach winding generated in the radial direction enters the 


shield material, 


mse round from 


EMG@w@tertal radial flux and iron Shield flux 


me- £ BtdcA=uy ff H.-aA (DS) 
S aA oO 4 O 
where 
— = 1 . — 7 = 
Hy (H aed ao 1. a (Ao eo thoaotlaag?) 16 (D6) 
a= rdddaz i. 
I ae as 1 o+ 4H aoe) H Seca: (D7) 
ao 2 rao xr sle(Or is) eS Riarcle, 2 “Gao 
fiemimmm1tS OL integration in the 9 direction will vary 
according to each element's included winding angle. Equation 
D5 may be written as 
De - Ros Puy ; Os paar eee eo merce OU (D8) 


Treating only the 


may be found from 


k=2, n=l, p=2 machine, the field strengths 


equation A&86 with the appropriate substitution. 


Once oe is determined, constrained aS mentioned above, 


the shield minimum thickness may be found as 


OR 


ies 








Figure D-2 Winding and Shield Configuration for Shield Thickness 
Determination 


ies 





Be Poe nooewrrrom lH Tgire loa, reference 9, to be 1.5 
Tesla. ‘The philosophy of shield design is based on the need 
tO protect personnel in an operating Space or prevent undue 


electromagnetic radiation from the vessel. 


Machine Losses 

Peenine losses may be attributed to two primary sources. 
Meeeees Such as bearing friction, windage, etc., are not 
addressed. 


meee conductor losses due to current flow, J. and J! 


Sone 
a C 2 2 pee Z 2p eee 
Po ax. OY Rag (ATR) FT SR 2") (D10) 





iMmeanation D110, the space winding factor, A, is assumed to be 


essentially the same for both armature and compensator 


windings. 
b) Eddy - current losses in the armature: 
poten. ° (=x) “oe 2d) oR = 
_ wa ao e ed’ w “av (D11) 
Peo a2 : 
where Bay is the mean-Saquare magnetic field seen by the 
ermMacure conductors: 
coe 
ly i 2 
Te : Bie 
Bay R -R_. J eee ( ) 
AO ai “RK 
ao 


Begs Given oy equations A56, A61, or AG3 as applicable. 


IS) 





Losses due to eddy currents and hysteresis losses in 
the shield is considered to be negligible because of 
the large armature-Shield separation (Ros <i Row)- 


Losses in the shaft material and armature winding 


Support structure are recognized but not addressed 


in this analysis. 
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Ae PN Dds 


DERIVATION OF MACHINE CHARACTERISTICS AND RATING EXPRESSIONS 


meeerring to Figure El, through trigonometric mani-ulation 


it may be shown that: 


Ep = Vin : x,t, isin Yo-— cos wvecot a) (El) 


The internal exitation voltage, Ep, May be expressed as: 


P 
DM ele 


5 a 


3) 
ll 

es 

NO 


Where, from equation B34 

















po po 
Pia e : Wa, _; wf 
fs = ly X_N AN, cl sin { 5 VS ert 5 ee ae 
NE 2 2 IR AO ae 
US ec eer tome ) (l-x*) ao 
(a3) 
and 
R 
a at + (Ee = 
ee ee (See oes 
C a Oo 2+p Ros 
jens (E4) 
p (4-p*) 
As a matter of convenience define 
ole pore 54. 
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FIGURE El. VOLTAGE-CURRENT RELATIONSHIP FOR 
OVEREXCITED COMPENSATED MACHINE. 
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Using equation £5, equation E3 may be reduced to: 

Mig = MIFA Noa ee (26) 
moe svOlt-ampere rating of the machine may be expressed as 
ol, = OE. fT, (=) (EW) 


From Appendix A, 


Nn, . modified 


(A 49) 


Substituting equations E2, E3, and A49 into E7, the machine 


rating may be expressed as: 











8 po We 

m4 7 Le ee b “4 Ewa ee gl 
Ee _- WU oJ I ¢8, (hoy Reo Cipestat 5) Sin (— on 
(E&) 


In order to solve for machine power from equation E8, an 


expression for (V7 ED) must be found. Define Ve and I, as 


rated machine base voltage and current (VEVip T,=1,)- Also, 


define the per-unit variables 
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< 
_ 
es 


ae. a FP 


ie Vi a Ip Vn 


Pevrading Cquation El by Ve ana the substitution of equations 


mo yields: 


Ng 


16 E Vi 





i (sin i ee VOXe ve ete sete!) (E10) 


With the per-unit synchronous reactance defined as 


SAB 


d Va 





Fquation E1O becomes: 


eve + x 


Cr fe Bt aes (Sima = =—=cos Ww Cot a) (E11) 


If equation Ell is divided by e, and simplified, the following 


a8 


relationship iS produced: 


= 


7 
Pi er 
il 
fea 
i 


x 1 (Sino cos cot a) Cee 


Where Xo the internally normalized per-unit synchronous 


feeactance is defined as: 








AB AA (212) 


Nee 





Semverting equation E12 from per-unit to absolute values 


yields: 


Ug 
a) ~ x_ (sin ~y - cos W cot a) (E14) 
Es 


Equation E13 may be expressed as: 


Circe ee a 0 
x, _ Sage ab’ “A (E15) 
os 


Using equation B25 for (Lo ~Loy)? and equations E2 and E3, 


equation E15 becomes: 











eee 
eee 3. sini) es C 
 -. a a @ | ao, 2p 1 sp (E16) 
eam Oe CUP wt. F6 (i-y*P) Sips 
Sin ( 5 ) 


All terms in equation E8 are now defined in terms of machine 
currents and geometry. Utlizing expressions developed in 
this and other appendices, it is possible to determine the 
compensator current density required for a given machine 
configuration and power level for any given value of a. 

Using equation C38 ,torque on the armature may be expressed 


as: 


2+k 2+k 











it = 24d 5 1 ) (80 ors ) ; Nees 
cee 2 2 en 
eooL 
KO ie KO 
[FL -Sin 5 )sin(ko) +7 2 sin ( 5 [eine (oa) | (C33) 
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where 2 and 2 are defined by equation A56. 
As a matter of convenience groups of terms in equation 


C38 may be defined as: 





eo a J Rao rc on 
7 | ee (E17) 
me O ' Ze AR 
Si 
NEe 
A. = J -2-sin (—z—) (Eee) 
iO 
A. = Ra Sai ( ) Coane) 
Equation C38 may now be written as: 
ae = A, [4,sin kero A,J sin k (d-a) | (E20) 
HereGwe on the armature may also be expressed as 
_ p (PWR) 
ta > OTF 
where (2) 


PWR = P cosy 


eiecemeuting equation E21] into E20 and simplifying, a set of 


@madradic equations in terms of sin (k®?) results: 


Q,sin* (ko) 4 Q,Sin (ko) + oe 0 een. 


BEY, 





Where, 


iy) on 2 mae | 3 

on =A) [do . A. +A. t20,A4,J cos (ka) ] (i230) 
p(PWR)A 

ee ap CN hg cos(ka) + 4.) (24 ) 


ep (pwR)* 


eae a Me? 
= eer, ays 
ae J 1 63 817 (ka) (EZoo 


Beege the equation for zero torque on the superconducting 


Pietra winding (C38 or C33) written as 
mesiko) = us A cos (ka) (E26) 


meeeon 22 reduces to: 





=e* \ a ce : = 
AgJ ah AAY AJ Aad + A, 0 (£27) 
where, 
nee ee l ‘ 
Ay Rag Res (ecucmtOms 4627) Or C30) 
A = (1923 } 
A, Le < Res (equations C34 or C35) 
a oe: ay 2 ee B29 
A. = Ay A, A“sin* (ka) ( ) 
Ne aa in’ (ka s (ka E30 
Ee = 2A A, A. A, sin“ (ka)cos (ka) ( ) 
p (PWR) A,A, W eoctko) sin (ke) 
— F2 2 2 in2z (ka) - ——$——$— $$ $$ 
A, = iy A, A, sin‘ (ka) TP 
_ ee oa: 
A, A, sin (ka) ( ) 


nets 





z A.A. A p(PWR) sin (ka) 
A. = - Z 
> = eee (32) 


Hes 


- p* (PWR)? 


A (E33) 
. An? PF? 
For a four pole machine with angle a fixed at 45° (90° 
electrically) with an iron shield boundary (assume Mie oss "009 tay 7m), 


first harmonic analysis of equations B29 through £33 yields: 


a ee A? 
a 0 
ee RA, *A° ~ ee (E34) 
eye 24, A. A (PWR) 
1 oe 
1) oh 


Beeeron £26 implies a set of conditions for zero torque on 
miemesuiperconauctinag field; one condition for each odd integer 
value of k. Higher order effects of the third and higher 
harmonics are not addressed herein. 


Consider equation E8 wirtten as 


V 
T 
JI J.(=) 
oa £ Le 


pa K (E35) 


IL sys. 





where, 





) 





0 pod 
_ 24 ee Be, uy a we wa 
Ky = (fe eo (l-y Reo yee | 5) Sin ( 5 


also, consider eacuation E16 written as 








ue 
cr — Ky J, 
where, 
2 sin( “way C 
2 ae ae —e 7 ; 5 a 
te oO ae. (yes) i Gewe 


Peeeeicuting equation E14 and E37 into E35 yields: 


J 


= = a 
pe= K Ji5, [1 Ky 7. 


l (sin Ww - cos Ww cot a)] 


faereton 639 may be reduced to: 


P= KJ - K 
a 


(sin wv - ont 
3 ae (Sin w COS. WcCOL. 0) 


4 


where 


ioe = KK 
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(E36) 


(B37) 


(B38) 


(Es 9") 


(F40) 


(E41) 





Using equations E7 and 40, an expression for J. in terms of 


J Y and ® may be found: 





i 
J - = Ke + KJ. (sin Ue IOS UWvCOt. oC) (E42) 
where, 
on eS) 
K,. = (E43) 
+) 
Tl ue 


If 4 is defined as sin J - cos J cot a, equation E42 may be 


[eae ten as: 


Neem (4 (E44) 





An examination of Figures 7, 8, and 9 reveals that the angle 
od between the armature phase a field and the superconducting 
main field is equivalent to the angle (~ + 90°) electrically 
at the zero torque condition. Therefore, equation E26 may be 


written as 
Sook ( + 90°) = J. A cos (ka) (E45) 
which reduces to 


Ceci) = a A cos (ka) (E46) 


141 


Angle ~ in equation E46 is a mechanical angular value 
Peabogous to angle ¢ in Figure 1. Employing the analytical 
Meemiiques which earlier resulted in equations E29 through E34, 
Mmemis alSO possible to obtain a fourth order equation in 
Sos sk © Or Sin k a as 


B,cos* (ka) + Byicos (scl) + B cos” (ka) + B.cos(ka) + B, = O 


3 is 


where Ba through Be are functions of machine geometry, Jos YW, 


and J. 
C 


Machine Length Parameters 
Based on machine lengths developed in reference 6, 


Table El represents the parameters used in this paper. 


iar 





Armature end turn length 


Armature self-inductance 
effective length 


Field-to-armature mutual 
inductance effective length 


Somoensator-to-armature 
mutual inductance effective 
feng th 


Eddy current loss effective 
length 


moncuction loss effective 
length 


TABLE El Effective Lengths 
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PROGRAM LISTING 


The program used in the development of this thesis, 
Meet Lf, accepts input data and provides output data as 
described in Table Fl. The program as listed herein provides 
output data for any given number of machine geometries and for 
each machine geometry provides data for armature current 
femetties ranging from .25 x 10° to 2.5 x 10° a/m*. Any 
input parameter may be incremented as required by insSerting 
meeroper @ Joop. The program accepts any number of pole 
pairs, any harmonic desired, and any value for shield 
permeability. Table F2 describes the arrangement of the 


input deck. 
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DATA ALGEBRAIC FORTRAN 


CLASS SYMBOL WORD 
Pnput p ie 
n N 
0 TWA 
wa 
0 TWC 
wc 
Of TWF 
a ALP 
1 F 
J JA 
a , 
Je JI 
af TA 
a 
HW . MUS 
cos wv PE 
ae RAI 
al 
R RAO 
aod 
Rey Gal 
se REO 
is cael 
om 
se RCO 
R_. Row 
Si 
R RSO 
so 
x L 
rar PE serene ws 7uey-tpyennmme nether womens 
Output all input quantities as above 
= C2 
1 B/G 
Cor C2is 
x XA 
a 
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ALGEBRAIC 
SYMBOL a 


DATA 
CLASS 


8 


ne 


ee 


Vt 


PWR 


Weert, Input/Output Data for THESIS 
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FORTRAN 
WORD 


xD 

(og 

NTA 

Wile 

APWR 

ee 

F4 

Pe 


Pek ieee wil 2. 


ae 


tee rogran 





DATA CARD # COLUMNS CONTAINS FORMAT 





ul 1-10 p Eom. 
11-20 N Jv Os: 
Zi 30 0 E103 
wa 
31-40 Q FilO.3 
WC 
A1-50 QO PO rs 
wet 
51-60 o} [SOG 8) 
2 1-10 J. lta 5 
11-20 J ¢ Be Otess 
5 1-10 Tn 1d IG Se 
A 1-10 i als 6: 
5 1-10 PF Pages 
6 1-5 IMUS LS 
6-10 KY iD 
y 1-10 Roi PaO. 3 
11-20 R PO ..3 
ao 
21.-30 Rea PaO 3 
31-40 Reo io. 3 
41-50 ae Plo. 3 
Gar 
51-60 R Pig. 3 
W 
61-70 Roi PLO. 3 
71-80 Re EOS 
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DATA CARD # COLUMNS CONTAINS POREATE 


8 a ee 0. I, EO 3 


TABLE F2. Arrangement of Input Data Deck 


Data cards 7 and 8 are sequentially repeated (without 
cards 1-6) depending upon the number of machines desired. 
KY is the number of machines being analyzed on each run. 
iemeiiay, Of course, be set artifically high and the program 
terminates once all data is processed. By setting IMUS = 
(eomecard 6, Oo 1S set equal to -l, assuming a laminated iron 
shield of high permeability; otherwise, IMUS = any interger 


not ¢ allows a to be computed. 
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